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ABSTRACT 

When considering the X-ray spectrum resulting from the reflection off the surface of 
accretion discs of AGN, it is necessary to account for the variation in reflected flux 
over the disc, i.e. the emissivity profile. This will depend on factors including the 
location and geometry of the X-ray source and the disc characteristics. We directly 
obtain the emissivity profile of the disc from the observed spectrum by considering 
the reflection component as the sum of contributions from successive radii in the disc 
and fitting to find the relative weightings of these components in a relativistically- 
broadened emission line. This method has successfully recovered known emissivity 
profiles from synthetic spectra and is applied to XMM-Newton spectra of the Narrow 
Line Seyfert 1 galaxy IH 0707-495. The data imply a twice-broken power law form 
of the emissivity law with a steep profile in the inner regions of the disc (index 7.8) 
and then a flat region between 5.6 Tg and 34.8 rg before tending to a constant index 
of 3.3 over the outer regions of the disc. The form of the observed emissivity profile is 
consistent with theoretical predictions, thus reinforcing the refiection interpretation. 

Key words: accretion discs - black hole physics - line: profiles - X-rays: general. 



1 INTRODUCTION 

In the canonical 'lamppost' model of AGN (George & Fabian 
1991), hard X-rays originate from a source in the hot corona 
surrounding the central black hole and are observed as a 
power law continuum in the spectrum of the object (the 
so-called power law component, or PLC). X-rays from this 
source will also be reflected off the accretion disc, with pho- 
tons being either Compton-scattered or undergoing photo- 
electric absorption which will lead to either Auger de- 
excitation or emission of a fluorescent line. This results in 
the so-called reflection-dominated component, or RDC of 
the spectrum. Within this component, a number of emis- 
sion lines from species in the disc will be observed (Fabian 
et al. 1989), as described by the self-consistent ionised re- 
flector model REFLIONX (Ross & Fabian 2005). 

The reflected flux from the accretion disc (the emisstv- 
ity profile) varies over the disc resulting from the variation 
in flux received from the coronal X-ray source as a func- 
tion of location on the disc. In the simplest case of a point 
source in flat, Euclidean spacetime, the emissivity at a point 
on the disc is proportional to the inverse-square of the dis- 
tance from the source, multiplied by the cosine of the angle 
at which the ray hits the disc from the normal, giving a form 

at large radius out from the source. In general relativity. 
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however, rays will be focussed towards the black hole (and 
the inner disc) so naively one would expect a steeper fall-off 
in emissivity with distance from the black hole. The details 
of the emissivity profile will depend upon the location, spa- 
tial extent and geometry of the source. 

When modelling the X-ray reflection spectrum from the 
accretion disc around the central black hole of an AGN, a 
(broken) power law emissivity profile (e cx r~°'), which de- 
scribes the reflected power per unit area from the disc, is typ- 
ically assumed and its parameters flt to the observed spec- 
tra. Power law forms of the emissivity proflle are motivated 
by calculations tracing rays from coronal X-ray sources in 
the spacetime around the central black hole on to the accre- 
tion disc (Miniutti et al. 2003; Dabrowski & Lasenby 2001; 
Suebsuwong et al. 2006). These calculations suggest a very 
steeply falling profile in the inner regions of the disc, then 
flattening off before tending to a constant power law slightly 
steeper than the law one would expect in the classical 
case. 

If the emissivity profile of an accretion disc can be de- 
termined from observations without a priori assumption of 
its form, this will allow for testing of the light-bending hy- 
pothesis and could place constraints on the location and 
geometry of the coronal hard X-ray source. 

We determine the emissivity profile of the accretion disc 
in the narrow line Seyfert 1 galaxy IH 0707-495, observed by 
Fabian et al. (2009). Previous attempts to recover the emis- 
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sivity profile of the accretion disc of MCG-6-30-15 from 
observed spectra by Dabrowski et al. (1997) fitted a free 
emissivity profile in the model of a relativistically broad- 
ened emission line, however here we recover the profile by 
considering directly the emission from each part of the accre- 
tion disc. The shape of a relativistically-broadened emission 
line, specifically the iron Kq line, is considered as the sum of 
contributions from successive radii in the disc and by finding 
the relative weightings of these components, the emissivity 
profile is recovered. 



2 RELATIVISTIC EMISSION LINE PROFILES 

From a given region of the accretion disc, the frequency of 
emission of an emission line as measured by an observer at 
infinity will be altered from the rest-frame frequency by rela- 
tivistic effects including Doppler shift, beaming and gravita- 
tional redshift, characterised by the transfer function (Laor 
1991), defined such that the spectrum of the line seen by the 
observer at infinity is given by 

Fo{uo)^ J Ie{re,^)T{re,g)dgredre (1) 

Where subscripts '0' denote quantities measured by the ob- 
server at infinity and 'e' those measured in the rest frame of 
the emitter. The redshift parameter is g — For a line at 
frequency i^e, the rest-frame spectrum is simply 

Ie{re,iy) = S{u - u^)e{r^) (2) 

The emissivity profile, e(r), gives the variation in reflected 
power as a function of radius, measured in the disc frame. 

It is clear from the form (1) that the line profile can 
be considered as the sum of contributions from successive 
radii in the disc, with the number of photons received by 
the observer at infinity from an annulus of the accretion 
disc radius Ve and thickness dre 

Noire, dre) = A'{re, dre)e{re) (3) 

Where A'{re,dre) is the projected area of the annulus as 
seen by the observer, which from (1) is the classical area 
of the annulus 2nredre, multiplied by the integral of the 
transfer function over all frequencies (or over redshifts, g, 
after a change of variables), where the transfer function is 
defined to give photon counts per unit frequency, rather than 
fiuxes (which is relevant here as when fitting to the reflection 
spectrum, the contribution from each region of the disc is 
given by its contribution to the total photon count). 



3 EMISSIVITY PROFILE DETERMINATION 

In order to determine the emissivity proflle of the accretion 
disc from the observed X-ray spectrum of an AGN, the rel- 
ativistically broadened iron Ka emission line (which has a 
rest frame energy of around 6.4 keV, though varies slightly 
with the ionisation state of the iron) is considered, as above, 
to be the sum of contributions from successive radii in the 
disc. A relativistically broadened line will have a different 
profile if emitted from annuli of different radius, which can 
be seen in Fig. 1, thus if the observed line can be broken 




Energy (keV) 



Figure 1. Relativistically broadened emission line profiles emit- 
ted from successive radii in the accretion disc ranging from 
1.235 rg to20rg and the integrated line over the disc out to 100 rg , 
calculated using the LAOR line model. 



down into contributions from different annuli and the rela- 
tive weightings of these found, the refiected fiux from each 
part of the disc (i.e. the emissivity profile) can be found. 

In the 3-10 keV region of the spectrum where the Ka 
line lies, the spectrum is well described as the sum of a 
power law (the PLC) and the refiected line (RDC) which 
is described by the reflionx model, once convolved with 
the KDBLUR model to apply the relativstic blurring in the 
accretion disc. The spectrum is modelled in xspec as 

powerlaw + ^ ^ kdblur ® reflionx (4) 

The radial intervals are selected to coincide with the radial 
binning of the transfer function used in the LAOR model 
(when fitting the model to the data, xspec looks up the 
transfer function from a pre-calculated table rather than 
evaluating it on-the-fly at each point, which would be com- 
putationally intensive) as it would be inappropriate to sub- 
divide the disc further, since the resulting line profiles from 
these subdivided bins would be essentially the same as the 
same transfer function is used. 

Only the 3-10 keV region of the spectrum is considered 
when determining the emissivity profile as this band con- 
tains a single emission line which can be decomposed into 
contributions identifying successive radii in the accretion 
disc according to the observed redshifts. When fitting for 
the emissivity, it is necessary to have a single emission line 
profile from which these components can be distinguished 
and as such the iron L region (0.5-1.0 keV), which contains 
multiple emission lines whose annular components will over- 
lap once relativistic blurring is applied, is not suitable as 
emission from different regions of the disc cannot be dis- 
tinguished here. Once the emissivity profile has been deter- 
mined, however, it can be applied to a reflection spectrum 
covering the full energy range. Furthermore, the EPIC pn 
detector on board XMM-Newton has a spectral resolution 
of 150 eV (FWHM) at 6 keV meaning there are 46 energy 
bins in the 3-10 keV band but between 0.1 and IkeV a res- 
olution of 100 eV (FWHM) results in only nine bins, with 
only three before emission lines from oxygen appear below 
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0.65 keV, thus in the iron L band, there are insufficient spec- 
tral bins available to fit for the emissivity profile. 

All model parameters are frozen to values which have 
been previously found to give the best fit for the source, ex- 
cept for the normalisations (photon counts) of the REFLIONX 
components, which are free parameters and are fit by min- 
imising to find the relative weightings of the components. 
The inner and outer radius of each annulus is set in the blur- 
ring kernel KDBLUR. Each annulus has constant emissivity 
(the index parameter is set to zero) such that there is no 
variation in emitted fiux across each annulus and the overall 
emissivity profile is given by the variation in photon counts 
between the annuli. 

The projected area of each annulus as seen by the ob- 
server at infinity is found from the transfer function in the 
LAOR model (which is used by kdblur for the convolution), 
by summing the photon count over all frequencies at a given 
radius with a constant, flat emissivity. The emissivity pro- 
file is then calculated by dividing the normalisations of the 
REFLIONX components by these areas as per Equation 3. 

3.1 Testing the Method 

It is not immediately obvious that a given emissivity pro- 
file will lead to a unique form of the observed spectrum, 
particularly once folded through the response function of a 
real X-ray detector, since the contributions to the line from 
different radii are not, strictly speaking, orthogonal func- 
tions. Before applying the method to real data, it is there- 
fore important to carry out preliminary tests, applying it to 
synthetic spectra whose emissivity profiles are known and 
checking that these profiles can be recovered correctly. 

This method has been applied to synthetic spectra in 
which line profiles with known emissivity profiles, formed 
from the LAOR emission line model, were folded through the 
response matrix of the EPIC pn detector on board XMM- 
Newton, using the fakeit command in XSPEC to simulate 
the spectra that would be observed if sources emitting such 
lines were observed with this instrument. Synthetic spectra 
were calculated for an exposure time of 300 ks and the inci- 
dent flux in the emission line was normalised to unity, since 
here we are interested only in the overall shape of the line 
spectra with no background. The response matrix of the de- 
tector maps only incident energies onto detector channels. It 
is independent of the incident fiux which will just provide an 
overall scaling that is irrelevant here since we are interested 
in whether the contributions from successive annuli can be 
recovered from the overall shape of the spectrum. 

A variety of line profiles were tested, with single power 
law emissivity profiles with indices ranging from 2.0 to 6.0. 
In addition to this, once-broken power law profiles were 
tested (these model line spectra are created using the LAOR2 
variant of the model). It can be seen from Fig. 2 that these 
profiles were successfully recovered. 

Scatter at larger radius appears where the radial step is 
finer than the binning used in the LAOR line model where the 
transfer function is evaluated. As a result, the line profiles 
from neighbouring radii in the same bin are essentially the 
same. They cannot be distinguished in the fit, leading to 
multiple points with essentially equal weighting rather than 
following the true emissivity profile as illustrated in Figure 
2(c). This scatter is reduced when the radial interval in the 
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Table 1. Values of frozen model pajrameters used in emssivity 
profile determination. 



Component 


Parameter 


Value 


powerlaw (PLC) 


Photon index, T 


3.09 




Photon count, norm 


1.50 X 10-3 


kdblur (RDC) 


Inclination, i 


53.96° 


reflionx (RDC) 


Photon index, F 


= powerlaw :r 




Iron abundance / solar 


8.88 




lonisation parameter, ^ 


53.44 




Redshift, z 


4.10 X 10-2 



1» 0.01 



E 10-= 



Energy (keV) 

Figure 3. Fit to the iron Ka edge (3-10 keV) of IH 0707-495 

with the disc reflection considered as the sum of components from 
successive radii. x^/DoF = 255.48/227 = 1.1255. 

fit is chosen to coincide with the radial binning of the LAOR 
transfer function. 



4 EMISSIVITY PROFILE OF IH 0707-495 

The emissivity profile of the accretion disc in this Nar- 
row Line Seyfert 1 galaxy was determined using the above 
method applied to the spectrum obtained using the EPIC pn 
detector on board XMM-Newton during four consecutive or- 
bits with an exposure time of 500 ks (Zoghbi et al. 2010). The 
values of the frozen parameters {i.e. those other than the 
photon counts from the reflection components) were taken 
from the best fit model of Zoghbi et al. (2010) and are shown 
in Table 1. 

The normalisations of the reflection components were 
initially fit in the energy range 3-10 keV, which encompasses 
the entire iron Ka line. The resulting fit is shown in Fig. 3 
and gives X^/DoF = 255.48/227 = 1.1255, indicating that 
the model describes the data well and that there are no 
significant components other than the power law continuum 
and refiection from the disc contributing to the spectrum in 
this energy range. 

The emissivity profile found from this fit by dividing 
the best-fit reflionx normalisations by the projected areas 
of the annuli is shown in Fig. 4(a). 

The model was also fit to the spectrum in the restricted 
energy range 3-5 keV. This energy range includes contribu- 
tions to the 6.4 keV emission line from regions of the disc 
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(a) (b) (c) 

Figure 2. Recovered emissivity profiles from synthetic spectra to test tlie metliod, each consisting of a single emission line with 
(a) a single power law emissivity profile, index 3 and (b) a once-broken power law emissivity profile, indices 6 and 2 with the break 
at a radius of 5rg. (c) is for the same profile as (b) but illustrates the scatter that occurs in regions where the radial points are 
more finely spaced than the radial bins in the LAOR model. 




a = 



1 10 100 1 10 100 

r/rg r/rg 

(a) 3-10 keV (b) 3-5 keV 

Figure 4. Emissivity profile of IH 0707-495 recovered from fitting to the spectrum in the energy ranges (a) 3-10 kcV and (b) 
3-5 keV. Error bars show la confidence limits corresponding to Ax^ = 1. Lines show power laws of index 6.0 and 3.3 in (a) and 7.8 
for the inner disc as well as a flat profile, for reference, through the middle region in (b). In the 3-5 keV fit, error bars are large at 
radii beyond 30 rg as emission in the iron Ka line from these radii is excluded from the fit range, so the points are unconstrained. 



out to 20 Tg , however contributions to the line emitted from 
radii beyond this are mostly excluded, with only their red- 
shifted wings included. Since the majority of the line flux 
is in the energy range above this and originates from the 
outer disc, the fit over the full line will be weighted heav- 
ily towards these radii, so fitting to the restricted 3-5 keV 
range will give better detail of the inner regions of the disc. 
The emissivity profile determined from the results of this 
restricted fit is shown in Fig. 4(b). 

From the above fits, a steep emissivity profile is ap- 
parent in the inner regions of the disc, with the 3-10 keV fit 



indicating a power law index of 6.0 out to a radius of around 
3rg and tending to a value of around 3.3 at large radius. 

In the 3-10 keV fit, the points in the region 4-20 rg are 
missing, with the best fit values to these normalisations be- 
ing zero. However the la error limits on these points are very 
large (on the arbitrary scale in Fig. 4(a), these points could 
be as high as 10~®), showing that the photon counts from 
this region of the disc are poorly constrained by these data. 
This is most likely due to the statistic in a fit over this en- 
ergy range being dominated by the line flux from the outer 
regions of the disc which is much greater than that received 
from the inner region. The line profiles from the innermost 
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Figure 5. Pseudo-image of the inner regions of the accretion flow 
in IH 0707-495 out to 60 rg formed from the observed emissivity 
profile. Shading corresponds to reflected flux in the disc frame on 
a logarithmic scale. 

regions out to 4rg account for the shape of the hne at low 
energy (and have a relatively high photon count due to the 
apparent steepness of the emissivity profile). Once these con- 
tributions and those from the outer disc are accounted for, 
the 4-20 rg region corresponding to the middle band of the 
emission line is not well constrained. This is remedied in the 
3-5 keV fit where the outer region of the disc (outside 20 rg) 
which dominated the line flux is removed and more detail of 
the inner disc is seen. 

The fit in the region 3-5 keV steepens the index of the 
emissivity profile to around 7.8 out to a radius of 5rg. This 
fit also indicates a flattening of the emissivity profile in the 
region 5-30 rg to an index of between and 1, however the 
error bars in this region are large and these points do not ap- 
pear in the 3-lOkeV fit because the contributions from these 
regions lie in the middle of the line profile so are masked 
by the contributions from either side. Note that in this fit, 
the points in the outer disc are poorly constrained and have 
large error bars because the energy range of their contribu- 
tions to the line is omitted from the fit. 

The determined emissivity profile can be used to con- 
struct a pseudo-image of the accretion fiow in IH 0707-495, 
showing the refiected fiux (measured in the rest frame of the 
disc) as a function of position in the disc (Figure 5). 

While this analysis was completed using the kdblur 
model which convolves the reflection spectrum with the 
LAOR line model for efficiency of computation, the analysis 
can be completed using other relativistic disc convolution 
models. One such model is the relconv model of Dauser 
et al. (2010) which convolves the refiection spectrum with 
the RELLINE profile, providing finer radial binning than the 
LAOR model. Using relconv yields very much the same 
result as is obtained here. A couple of further points are 
obtained in the flat 5-30 rg region of the emissivity profile 
when performing the fit over the full 3-10 keV range consis- 
tent with the above result, however this region is still largely 
unconstrained (until the 5-10 keV band is excluded from the 
fit). This indicates that this is not due to systematics in the 
LAOR model but is rather the masking of this region by the 
high flux from the outer disc in the observed spectrum. 



Table 2. Fit to the emissivity proflle of IH 0707-495 (3-10 keV) 
using the kdblurS model (twice- broken power law emissivity). 



Parameter 


Fit Range 


Value 


Error (Icr) 


Index 1 


2-12 


7.83 


-0.66,-1-3.97 


Break radius 1 


1.235 - 10 rg 


5.60 rg 


-0.48,-1-0.15 


Index 2 


0-2 


7.84 X 10-5 


-1-0.38 


Break radius 2 


5 - 400 rg 


34.75 rg 


-4.72, -1-4.42 


Index 3 


1.5 - 10 


3.30 


-0.32, -1-0.43 



5 TESTING THE RESULT 

It is important to test these results for the emissivity proflle, 
in particular the flattened region of the proflle at 5-30 rg, by 
fitting one continuous refiection spectrum with the appro- 
priate emissivity profile over the entire disc. It would appear 
from the above that the closest analytic form for the profile 
is a twice-broken power law. 

In XSPEC, the emissivity proflle of the disc is applied 
by the relativistic blurring kernel, kdblur (single power 
law emissivity) or kdblur2 (once-broken power law) , to the 
rest-frame reflection spectrum reflionx. The blurred spec- 
trum is obtained by convolving the original spectrum with 
the proflle of a relativistically broadened line (delta-function 
spectrum in the rest frame) from the LAOR or LAOR2 models. 

Modifled versions of these models were written to give 
a relativistically blurred reflection spectrum from the disc 
corresponding to a twice-broken power law emissivity pro- 
file: The model kdblur3, which provides the blurring kernel 
based upon the line profile laor3^, with which xspec may 
convolve the rest frame spectrum. These models have the 
three power law indices and two break radii for the emissiv- 
ity profile as input parameters. 

The spectrum was fit in the region 3-10 keV (the iron 
Kq line, with no significant spectral components other than 
the power law contirmum and reflection component) with 
the model 

powerlaw + kdblurS ® reflionx (5) 

with the power law indices and break radii of the emissivity 
proflle flt as free parameters, within limits to constrain them 
to roughly the values observed (i.e. a steeper index tending 
to a moderate index, with a flat proflle in the middle region) 
but sufficiently broad not to place any particular bias on the 
results. Other parameters of the models were again set to the 
best-fit values of Zoghbi et al. (2010). The results of this fit 
are shown in Table 2. 

Table 3 shows this fit compared to an equivalent fit 
using only a once-broken power law as had been assumed in 
previous work in the analysis of the spectra. 

Modelling the emissivity profile as a twice-broken power 
law fits the data substantially better than using a once- 
broken power law, giving a reduction in oi 23.2 (a re- 
duction from 1.11 per degree of freedom to 1.04 per degree 
of freedom). The best-flt parameters obtained using the kd- 
blurS model agree well with the emissivity proflles recov- 
ered from the sum of reflection components, with a power 
law index of 7.8 in the inner disc out to 5rg and flattening 

^ http:/ /www-xray.ast. cam. ac.uk/-drw/models/kdblur3 
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T il l <-) 

kablur2 


kdblurS 


Index 1 


4.82 


7.83 


Break radius 1 


6.85 rg 


5.60 rg 


Index 2 


2.09 


7.84 X 10-5 


BreaJc radius 2 




34.75 rg 


Index 3 




3.30 




287.52 


264.32 


xVDoF 


1.1144 


1.0365 



Table 3. Comparison between fits to the iron Ka line (3-10 keV) 
of IH 0707-495 using models with once- and twice-broken power 
law emissivity profiles (kdblur2 and kdblur3, respectively). 



0.1 - 

0.01 - 

0.001 - 

0.0001 - 

1e-06 - 

1e-06 - 

1e-07 - 

1e-08 - 

1e-09 - 

1e-10 — 
1 



kdblura 
kdblur2 



Figure 6. Comparison between the best fit emissivity profiles 
obtained for IH 0707-495 using the KDBLUR3 (x^/DoF = 1.037) 
and KDBLUR2 (x^/DoF = 1.114) models. 



to a very low index of 8 x 10"^ from here to 35 rg before 
tending to an index of 3.3 over the outer regions of the disc. 
The best fit emissivity profiles obtained using the KDBLUR3 
and KDBLUR2 models arc shown in Fig. 6. 

The parameters are well constrained to these values by 
the 1(7 confidence limits, with the exception of the index in 
the middle, flattened region, where the lower limit is not well 
constrained by the data but an upper limit of 0.4 is placed 
on the index here. Likewise, the lower limit of the power law 
index for the inner disc is found, yet the upper bound is less 
well constrained with a confidence limit of +3.97. This is 
likely due to the low flux received from the inner disc since 
photons originating from iicrc arc more likely to fall into the 
black hole, thus increasing the emission from the inner disc 
by steepening the emissivity profile will have less of an effect 
on the line profile detected at infinity. 

The fits to the spectrum using the two models arc shown 
in Fig. 7. The form of the edge between 6.5 and 7keV is 
reproduced much better by the KDBLUR3 model than by 
KDBLUR2. 



5.1 Further Spectral Components 

Adding a narrow (Gaussian) line at 6.4 keV (fitting for the 
line width, a, and photon count) to account for X-ray emis- 
sion from a more distant refiector from which relativistic 
blurring is not seen does not significantly improve the fit to 



the data. The reduced changes to only 1.01 (compared 
with 1.03) when the emissivity profile is frozen at the val- 
ues determined above. The fiux in the narrow line is 100 
times less than that seen in the broad line from the disc, 
indicating that this component is not significant. Allowing 
the parameters of the kblur3 emissivity profile to vary also 
in the fit yields a similar result with the observed emissivity 
profile not changing greatly (it is steepend slightly over the 
outer disc to compensate for the extra flux in the narrow 
line coinciding with energies from the outer disc, but the fit 
is not improved with a reduced of 1.03). 

Adding a second ionisation component in the disc (ei- 
ther simply adding a second component that is cospatial 
with the first, or dividing the disc into two ionisation com- 
ponents for the inner and outer disc and fitting for the break 
point) does not improve the fit. The reduced x^ statistic is 
unchanged at 1.03 when the disc from 45 rg outwards has 
the ionisation parameter reduced to ^ = 28 (while that in 
the inner disc is still found to be around ^ = 53). The best 
fit emissivity profile for the disc with two ionisation compo- 
nents is found to be essentially unchanged from the previous 
result using one component. 



6 DISCUSSION 

It is evident from the results here that the emissivity pro- 
file of IH 0707-495 resembles the form predicted by Miniutti 
et al. (2003), with a steep power law (e oc r"™) in the in- 
ner regions of the disc, found to fiave an index of around 
7.8, fiattening to a very shallow index, close to zero, before 
tending to a constant index of 3.3, slightly steeper than the 
Euclidean case, at large radius. 

Such an emissivity profile suggests that 90% of the X- 
ray flux reflected from the accretion disc (where F{< r) = 
f e(r') r'dr', measured in the local disc frame) is reflected 
from the innermost regions of the disc, within 2rg of the 
central black hole as shown in Fig. 8 (the solid line shows 
the cumulative flux distribution measured in the disc frame 
and the dashed line that measured by an observer at infinity 
where the outer disc is more dominant due to the larger solid 
angle subtended and the tendency of photons emitted from 
the inner regions of the disc to fall into the black hole — 
£is seen at infinity, only around 66% of the refiected fiux 
originates from the inner disc, within 4rg). 

Qualitatively, this form of the emissivity profile can be 
understood by considering an isotropic point source above 
the disc in a fiat, Euclidean spacetime. The fiux received 
from the source at each point on the disc, and thus the 
reflected power (emissivity) from that point will vary simply 
as the inverse square of the distance from the primary X-ray 
source, projected into the direction normal to the disc plane, 
i.e. if the source is at a height h above r = in the disc, the 
emissivity will go as (r^ + h?)^^ cos d. with cos d = , 

the angle from the normal at which the ray hits the disc. 
This will be constant (a flat profile) in the limit r <^ h (the 
inner disc) and tending to r~^ at large radius, giving a power 
law emissivity profile with index 3. 

In the presence of the black hole, gravitational light- 
bending will act to focus light rays from the point source 
towards the black hole and will increase the flux incident 
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Figure 7. Best fits to the iron Ka edge of IH 0707-495 using (a) twice-broken and (b) once-broken power law emissivity profiles 
(kdblurS and KDBLUR2 respectively). 



Disc frame - 
Observer at infinity 



Figure 8. Cumulative refiected X-ray flux distribution, F{< r), 
as a function of radius in the disc as implied by the determined 
emissivity profile for IH 0707-495 measured in the disc frame com- 
pared with that as measured by an observer at infinity (taking 
into account the transport of radiation from the disc to the ob- 
server around the black hole for IH 0707-495 which is observed at 
an inclination angle of 53°). 



on the inner regions of the disc relative to the outer re- 
gions. This will cause the emissivity to fall off more steeply 
in the central regions of the accretion disc as is observed with 
IH 0707-495, before reaching a power law index of around 
3.3 at large radius, slightly steepened from the classical value 
of of 3.0. 

In addition to gravitational light bonding focussing the 
light rays on to the disc, the emissivity profile will be affected 
by relativistic effects on the disc itself. While the classical 
area of an annulus radius r and thickness dr is 2'nrdr, the 
proper area as measured in the rest frame of the disc for a 
given dr is increased as r decreases, as the space is warped 
close to the black hole and proper radial distances increase 
at small r as described by the Kerr metric. As the disc or- 
bits the black hole (which on the last stable orbit will be 
at half the speed of light) , the disc material will be length- 
contracted by a factor of 7, the Lorentz factor according to a 
stationary observer, so the area of the radial bins is increased 



further in the disc frame by a factor of 7. These eflects in- 
crease the proper area of a radial bin closer to the black hole 
so will tend to reduce the steepness of the emissivity profile 
(as the incident flnx is divided by the disc area). 

The dominant effect steepening the emissivity profile 
close to the black hole is dilation of the proper time in the 
disc frame (i.e. gravitational redshift). If the X-ray source 
emits photons at a constant rate in its own rest frame, the 
arrival rate in the disc frame will be enhanced if that disc 
element is closer to the black hole, as its proper time elapses 
more slowly. To an approximation, the photon arrival rate 
varies between different regions of the disc like t{r), the 
derivative of co-ordinate time with respect to the disc el- 
ement's proper time, which increases steeply moving closer 
to the black hole. Furthermore, since emissivity is defined 
as the flux emitted, a further factor of approximately i is 
introduced as each photon is blueshifted as it travels toward 
the black hole. 

If the X-ray source is moving {e.g. orbiting the rota- 
tion axis), emission will be relativistically 'beamed' into the 
direction of motion of the source. This will enhance emis- 
sion in front of the source while decreasing omission behind 
and to the sides of the direction of motion. If the source is 
orbiting the rotation axis on a locus above the disc plane, 
emission will be enhanced on to regions of the disc below 
the orbital locus, while the emissivity outside the orbit is 
reduced. If the source is close to the black hole, this will 
further steepen the emissivity profile. 

The combination of these effects is illustrated by theo- 
retical emissivity profiles created by a code written to trace 
rays from an isotropic point source in the corona above the 
plane of the accretion disc and counting the photons inci- 
dent on radial bins in the equatorial plane, dividing by the 
area of the bins while taking into account the above effects 
(Fig. 9). Such simple simulations can qualitatively explain 
the form of the observed emissivity profile. 

The best fit to the spectrum with reflection components 
originating from successive radii indicates a non-zero contri- 
bution to the spectrum from photons originating from the 
inner regions of the disc right down to the innermost stable 
orbit for a maximally spinning Kerr black hole at 1.235 rg. 
This suggests that the accretion disc extends down to the 
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Figure 9. Theoretical accretion disc emissivity profiles due to (a, 
black) a stationary point source on the rotation axis above the 
disc plane, (b, blue) and (c, red) ring sources co-rotating with 
the disc below, from ray tracing simulations, taking into account 
relativistic effects on the area of the accretion disc. The emissivity 
profile steepens over the innermost regions of the disc to a power 
law index of up to 7 in (a) and 5 in (b) then flattens off before 
tending to a constant power law index of around 3 over the outer 
disc. Increasing the radius of the ring source moves the break- 
point between the flat region and the outer disc outwards to the 
point approximately below the source, (c) illustrates the extreme 
case of a source at large radius while at a low height above the 
disc, enhancing the steepening over the inner regions of the disc. 
The emissivity here rises directly below the source with a steep 
fall-off at large radius since due to the low height, many photons 
are intercepted by the disc before they can propagate further. 



innermost stable orbit. The contribution received from the 
innermost regions of the accretion disc by an observer at 
infinity is, however, small once the projected area of the an- 
nuli (including the effect of the propagation of rays around 
the black hole to the observer) is accounted for. As such, the 
drop in the observed emissivity over several decades across 
the innermost regions may not be as significant to the ob- 
served line profile as first it seems. 

The steep power law emissivity profile determined here 
without a priori assumption of its form and the agreement 
of this form with theoretical emissivity profiles calculated 
for point sources around the black hole show that the steep 
emissivity profiles previously required to explain the emis- 
sion seen in terms of reflection of X-rays from a hard X-ray 
source near the black hole off of the accretion disc are indeed 
plausible. The refiection spectrum has a significant contri- 
bution from the inner regions of the disc right the way down 
to the inner-most stable orbit. 

While the profiles of emission lines originating from 
annuli in the accretion disc are calculated using the LAOR 
model and thus implicitly assume a maximally spinning Kerr 
metric (a — 0.998), the profile of the line from a general 
annulus does not vary significantly with the spin parame- 
ter. This can be demonstrated by computing line profiles 
outwards of 6rg (i.e. the innermost stable orbit of a non- 
spinning, Schwarzschild black hole) for varying values of the 
spin parameter using a relativistic emission line model such 
as RELLINE (Dauser et al. 2010). The effect of varying the 



spin parameter is to move the innermost stable orbit, defin- 
ing the innermost possible extent of the accretion disc and 
thus the extremal gravitational redshift, as illustrated in Fig- 
ure 10. As each annulus is considered independently, varying 
the black hole spin will determine whether contributions are 
detected from the inner annuli whose existence depends on 
the location of the ISCO. As such, the spin of the black hole 
will not affect the determined emissivity profile for the outer 
disc and the inference of emission from regions of the disc 
down to 1.235 rg implies a spin of the black hole close to 
maximal. 

Theoretical emissivity profiles calculated here are for 
the simple case, considering only direct radiation from a 
point source. We have not in the present work considered the 
effect of spatially extended X-ray sources or the effect of ra- 
diation returning to the disc after reflection and hence being 
reflected for a second time (Cunningham 1976) which may 
modify the final emissivity profile (for instance, increasing 
the reflection from the inner regions of the disc as more rays 
will be focussed to regions closer to the black hole, giving a 
steeper inner emissivity profile) and sharpen the spectrum 
(Ross et al. 2002). This will be considered in future work. 

It should be noted that if this twice-broken power law 
for the emissivity profile is averaged over the disc to a single- 
slope power law (averaging the index over radius, weighted 
by the annular area, rdr), the single power law index is 
found to be around 3.3. This would be found if a refiection 
component blurred using a single-slope emissivity profile is 
fit to the spectra and is only slightly steeper than the clas- 
sical case, masking the significant steepening over the inner 
disc by relativistic effects. 

Throughout the preceding analysis, the standard limb- 
darkening of emission from the accretion disc (Laor 1991) is 
assumed. Work by Svoboda et al. (2010) suggests that the 
disc emission due to refiection can in fact be limb-brightened 
or at least isotropic from the innermost regions of the disc 
close to the central black hole. This will lead to an enhanced 
red wing of the emission line over the typical limb-darkened 
profile as can be seen in the relline model of Dauser et al. 
(2010), and as such will lead to an apparently steeper emis- 
sivity profile. Thus, in reality the power law slope over the 
inner disc may be less than 7.9 in IH 0707-495, while the ef- 
fect on the emissivity profile obtained will be less relevant for 
the outer disc regions where relativistic effects are weaker. 



7 CONCLUSIONS 

The X-ray reflection emissivity profile of the accretion disc 
in IH 0707-495 has been determined from the proflle of the 
iron Ka emission line in the X-ray spectrum without a prion 
assumption of its form, by considering independently the 
relative contributions from successive radii in the disc. 

The emissivity profile was found to agree with theoret- 
ical predictions; a steeply falling profile in the inner regions 
of the disc with a power law index of 7.8, then flattening to 
constant emissivity between 5 rg and 35 rg before tending to 
a constant power law of index 3.3 over the outer regions of 
the disc. The profile obtained suggests emission right down 
to the last stable orbit at 1.235 rg. 

The emissivity profile obtained is consistent with fit- 
ting a relativistically blurred reflection spectrum from a disc 
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with a continuous twice-broken power law emissivity pro- 
file, which provides a better fit than the previously assumed 
once-broken power law. 

The observed profile can be understood in the context of 
gravitational light bending as rays travel around the black 
hole from the source to the disc before being reflected, as 
well as relativistic effects on the disc and agrees with forms 
predicted by ray tracing simulations, indicating that reflec- 
tion from an accretion disc with a steep central emissivity 
profile is plausible. 

The results can be compared to theoretical emissivity 
profiles computed for a range of locations and geometries of 
the X-ray source as well as disc geometries in order to place 
constraints on the properties of the coronal X-ray sources in 
AGN. 
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Figure 10. Theoretical line profiles for varying values of the black 
hole spin parameter, a, calculated using the relline model, (a) 
shows emission lines emitted from a region of the disc extending 
from the relevant innermost stable orbit (ISCO) for the value of 
a out to 400 rg, with the extended redshifted wing of the line 
for higher values of a resulting from the accretion disc (and thus 
emission) extending down to an ISCO closer to the black hole, 
(b) and (c), showing lines emitted from the regions of the disc 10- 
400 rg and 10-20 rg respectively, demonstrate that the effect of 
the spin on the emission from a given radius is small. The effect 
of spin decreases in significance further from the black hole. 
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